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Abstract 

It is stated that holonomy corrections in loop quantum cosmology intro- 
duce a modification in Friedmann's equation which prevent the big rip singu- 
larity. However in [1 1 it has been proved that this modified Friedmann equa- 
tion is obtained in an inconsistent way, what means that the results deduced 
from it, in particular the big rip avoidance, are not justified. The problem is 
that holonomy corrections modify the gravitational part of the Hamiltonian of 
the system leading, after Legendre's transformation, to a non covariant Lag- 
rangian which is in contradiction with one of the main principles of General 
Relativity. A more consistent way to deal with the big rip singularity avoid- 
ance is to disregard modification in the gravitational part of the Hamiltonian, 
and only consider inverse volume effects |2 1. In this case we will see that, not 
like the big bang singularity, the big rip singularity survives in loop quantum 
cosmology. 



Pacs numbers: 98.80.Qc, 04.20.Dw, 04.60.Pp, 04.20.Fy 

1.- Introduction — Observations from the spacecraft Wilkinson Microwave Aniso- 
tropy Probe (WMAP) indicate that our universe would be dominated by "phantom 
energy". This fact would be explained assuming the presence of dark energy with 
a negative pressure, which could be modeled by a barotropic fluid with state equa- 
tion P = Lop, with uj = —1.10 lb 0.14 (see [3]). However, a universe containing 
dark energy modeled by equation P = ujp with a; < — 1 lead to a classical future 
singularity called big rip ||4]|. 

This singularity is at classical level, i.e., it appears solving the classical Raychaudury 
equation. But, if one took into account quantum effects it could be avoided, that 
is what happens in loop quantum cosmology where a good number of papers have 
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shown that this singularity is replaced by a non-singular bounce (see for instance 

mmiziiiiii). 

The main idea is that loop quantum cosmology assumes a discrete nature of 
space, which leads, at quantum level, to consider a Hilbert space where quantum 
states are represented by almost periodic functions of the dynamical part of the 
connection. Unfortunately, the connection variable doesn't correspond to a well 
defined quantum operator in this Hilbert space, what leads to re-express the grav- 
itational part of the Hamiltonian in terms of almost periodic function, which could 
be done from a process of regularization. This new regularized Hamiltonian, which 
contains the quantum geometric effects, introduce a quadratic modification (p^) in 
the Friedmann equation at high energies which give rise to a bounce when the en- 
ergy density becomes equal to a critical value of the order of the Planck energy 
density. 

Recently, we have shown in |[T]| that the regularized Hamiltonian, after Le- 
gendre's transformation, gives rise to a non-covariant Lagrangian, which contra- 
dicts the main principles of General Relativity. This means that the modified Fried- 
mann, which is obtained from this regularized Hamiltonian, is incorrect. Thus, the 
current statement that, in loop quantum cosmology, the big rip singularity is re- 
placed by a non-singular bounce, is wrong. 

To sort out this situation, it seem clair that one cannot change the gravitational 
part of the Hamiltonian (or it has to be changed by another one that doesn't lead 
to a non covariant Lagrangian), which means to disregard the Hilbert space of al- 
most periodic functions. In this way we review the early Bojowald's work where 
the author only takes into account inverse volume correction in the energy dens- 
ity. These corrections lead to a new modified Friedmann equation where the big 
bang singularity is avoided, but the big rip one survives. Unfortunately, to main- 
tain, at the same time, the discrete nature of the space and the gravitational part of 
the Hamiltonian it seem that one has to choose a Hilbert space where the Hubble 
parameter, or the connection, takes values on a circle rather in the real line, which 
seems, from a physical viewpoint, an unjustified assumption. 

The units used in the paper are h = c = Mp = 1 being Mp = (g^)^^^ the 
reduced Planck mass. 

2. - Classical cosmology and the big rip singularity — For the flat FRW space-time 
filled by a perfect fluid with state equation P = up (being lu < —1), Einstein 
theory is obtained from the Lagrangian C = \Ra^ - pa^ where R = 6{H + 2H'^) 
is the scalar curvature, a is the scale factor and H = ^is the Hubble parameter. 

This Lagrangian can be written as follows £ = 3 ^ "'^^^ — a?aj — pa?, this 
means that the same theory is obtained avoiding the total derivative, which gives 
the Lagrangian Ce = —iH^a? — pa^. The conjugate momentum is then given by 
p = = —6Ha^, and thus the Hamiltonian is 

-He = ap- Ce = -SH'^a^ + pa^- (1) 
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In general relativity the Hamiltonian is constrained to be zero which gives the 
Friedmann equation 



= p/3. (2) 
This equation with the conservation equation 

d{pa^) = -Pd{a^) ^ p= -3H{1 + u;)p, (3) 

are the dynamical equations that describe the evolution of our universe. 

Taking the derivative of the Friedmann equation with respect to the time and 
using the conservation equation one gets the so-called Raychaudury equation 

H = -^{l + cv)H^ (4) 

that together with the Friedmann equation give the solution 

2 1 4 1 

= 3{l+u)t-trip 3(l + w)2(i-is)2' 

where trip = to — 3Ho{i+uj) '•^^ which the big rip is reached, and Hq = 

H{tQ) is the current value of the Hubble parameter. 

Since observations from WMAP give the approximate values Hq = 2.3 x 
^^^^^ seamd ^^'^ ^ ~ ^^^^ mcans that the time from the present to the big 

rip is approximately 90 bilhon years. 



3. - Loop quantum cosmology — The old quantization of loop quantum cosmology 
was done using two canonically conjugate variables, one of them was the dynam- 
ical part of the connection, namely c, and the other one was the dynamical part 
of the triad, namely p. These variables are related with the scalar factor and the 
extrinsic curvature K = ^ahy the relations (see for instance ifTOl ) 

IpI = a^, c = 2jK = 'ja, (6) 

where 7 is the Barbero-Immirzi parameter, and their Poisson bracket is given by 

{^>P} = i|-|l^ = i^gn(p)- 

Alternatively, one can introduce the new canonically conjugate variables ifTTl 

^ = ^ = 1^1 = IPl'^' = ("7) 

with Poisson bracket V} = ^sgn(V). 

To built the quantum theory, in LQC, is usual to choice as a Hilbert space, the 
quotient space to the Besicovitch space of almost periodic functions by its subspace 
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of null functions. The Besicovitch space (see for details |fT3l ) is the closure of 
trigonometric polynomials under the semi-norm (in the c-representation) 

11*11'= lim ^ r |*(c)|'dc, 
where c is the connection. And an element of this space as the expansion 

with //„ G M and a„ G (the space of square-summable sequences). 

Remark 0.1. It is important to realize that ||.|| is a semi-norm. For example, all 
continuous functions with compact support have zero norm. 

In this space on can define the operator p as follows p = — However 
the operator c defined by c^'(c) = c^'(c) is not well defined in this Hilbert space 
because for a general quantum state ^'(c) = Yln&'L ctnC*^"'^''^ one has 

Nc^lP = lim — > lonP = +00. 
L-s>oo 3 ^ 

Note that, one can also use the /3 representation, where now the semi-norm is 
||*||2= lim ^ r 
and the expansion of quantum states is 

^{fi) = J2 O^nWn) = Yl ^ne""-^^^. (8) 

In this representation the operator V is given by —YZp' ^'^^ same 

reason as c, is not well-defined. 

Now the volume operator, namely V, is given by = | V|, and act as follows: 

Remark 0.2. It is common, in loop quantum cosmology, to use the v-representation 
that could be defined as follows: 

Then, for a general state defined by dS| one has '^{v) = Cin^{v — z/„). 
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The operator V, in this representation, acts as follows: 

where we have used integration by parts. And thus, V'^{v) = 
Finally, the operator e^"l^ acts as follows: 

^^(i/) = — [ e'''^^{l3)e-'''^l'^dl3 = ^{v - 2a). 

We can also define negative powers of V in the following way: First we use 
the formula 

A7 

where A is a parameter with dimensions of length, which is determined invoking the 
quantum nature of the geometry, that is, identifying its square with the minimum 

eigenvalue of the area operator in LQG, which gives as a result A = y^'^T (^^^ 

ini). 

Then, from this formula, we define the operator |V|^^/^, which applied to the 
eigenstate gives 

= ^sgn(V) (f^[e^, |V|^/2] ^ ^^^^ |V|i/2]e;^) \u) 



-^\\u + 2X\^/^ -\u-2X\^/^\\i^). (9) 
AV7 



Note that for large values of one has 

1 111, + 2A|i/2 2A|i/2|| - —W\-'/^ 

which is the inverse of the square-root of the eigenvalue of the volume operator. 
And for small values of v on has 

1 V + 2A|l/2_|^_2;^|l/2p, 1 



A^ V27A3/2' 

which vanishes a.tv = 0. 

Remark 0.3. Note that there is an ambiguity in this procedure because one could 
have started with \V\~^ where j G (0, 1) in place oflVl"^^"^ (see [11]). 

Here it appears a very important problem, because in the gravitational part of 
the Hamiltonian it appears /? or c. In fact, one has 
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and since the operators fi and c are not well-defined, in order to built the quantum 

theory, one needs to re-define the Hamiltonian. To be precise, we will work in 

_ -M 

the P representation and first at all we consider the holonomies hj{\) = e * 2 "^j, 
where aj are the Pauli matrices. 

Then, since doesn't have a well-defined quantum operator, to construct a 
consistent quantum Hamiltonian operator, one needs one almost periodic function 
that approaches for small values of /3. This can be done using the general 
formulae of loop quantum gravity to obtain the regularized Hamiltonian 



xhj\X)hk{X){h-,HX),\V\} 



+ p\V\, (10) 



which captures the underlying loop quantum dynamics (see ifTOl [T4l ). Note that, 
this Hamiltonian admits a quantum version because the operators hj{X) are well 
defined in the Hilbert space. 



4. - Effective dynamics — Coming back to the regularized Hamiltonian ( 10 1, an easy 
calculation gives |[T5l IT6II 

^LQC = -3|V|^^^^ + /,|V|. (11) 

Now, the Hamiltonian constraint is then given by ''^'^^•^2^^ = |, and the Hamilto- 
nian equation gives the following identity: 

V = {v,nLQc} = "i^^^gnW ^ (12) 

H = ^ /5 = ^ arcsin(2A7F). (13) 



Writing this last equation as follows H = (1 — sin (A/3)) and using 

the Hamiltonian constraint H lqq = <;=^ '''"/A^'* = | one obtains the following 



modified Friedmann equation in loop quantum cosmology 

being pc = This equation with the conservation equation p = —3H{1 + u})p 
gives the dynamics of the universe in loop quantum cosmology. 

Integrating these equations one can obtain that if the current energy density of 
our universe is p{to) = po, then the dynamics of the universe will be given by 



p{t)=[-{l + u;yit-ty + -] , H{t) = ^^it-t)pit) (15) 
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where t = to 7^ — T^^^^ solution is defined for all time and it finishes 



at {H = 0,p = 0). The main property of this solution is that the universe is 
expanding phase until time t = i. At this time it bounces and re-collapses forever 
and ever. 



It is clear that the non-singular behavior described by equation ( 15 I, is very 
different to the classical one where the universe presents a big rip singularity. 

However, as we have stated above, loop quantum cosmology could be built us- 
ing two canonically conjugate variables, (V, /3) which are related with the standard 
variables through the relations |V| = and f3 = -jH (see for instance formulas 
2.1 and 2.2 of |17|). However, if one starts directly from the effective Hamilto- 
nian (111, although it is assumed that |V| = a^, the definition of the variable /3 is 



obtained from the equation V = {V, 7^lqc}> which gives = jx 8i,Tcsm{2XjH) 
and differs from the initial definition of (3. Consequently, if one takes |V| = 
and P = jx ai'csin(2A7i/) as a canonically conjugate variables, then the standard 
variables {a,p = —QHa?) do not remain canonically conjugate, because now 

= x/r^^4^^A27P / constant. (16) 

Moreover, it's well-known that the other current cosmological theories are built 
from two invariant, the scalar curvature i? = 6 (^H + 2H^^ and the Gauss-Bonnet 

curvature invariant G = 24:H^ (^H + H^^ (note that the Weyl tensor vanishes 
for FRW cosmologies). For example, in Gauss-Bonnet gravity lITSl it is used the 
Lagrangian Cmg = o,'^f{R, G) — a?p, and semiclassical gravity, when one takes 
into account the quantum effects due a massless conformally coupled field (see 
for instance L19J ). is based in the trace anomaly T^ac = aDi? — f G (being a > 
and /3 < two re-normalization coefficients). However, from the Legendre 
transformation Hlqc = — ^sgn(V)V/3 — Clqc one gets, in terms of the standard 
variables, the following Lagrangian 

Clqc = avcsm{2XjH) + ^ (l - y^l - Aj^X^^j - a'p, (17) 

which is not invariant. This is in disagreement with one of the main principles of 
general relativity. 



From these observations, one can conclude that using the Hamiltonian (111, 
some contradictions like the variables (a,p) are not canonically conjugates or that 
the Lagrangian is not invariant will appear. In particular, the modified Friedmann 



equation was obtained assuming that (111 is the Hamiltonian of the system, and 



thus, one can conclude that this equation has been obtained in an incorrect way. 

The question is, why one uses this Hamiltonian in LQC? The answer seems 
clear: since the operator is not defined in the Besicovitch space of almost peri- 
odic functions on has the replace the standard Hamiltonian by another built with 
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variables that admit a well-defined quantum operator in this space. However, an- 
other viewpoint is possible, one can use a different Hilbert space where was 
well-defined. One of these space was introduced by Bojowald in ||20ll2T1l . it was 
C'^{SU (2), diiu) the space of functions of isotropic connections which are square 
integrable with respect to Haar mesure. Essentially, this space, in the c representa- 
tion, is the clause of continuous functions with the norm 



^r||2 = / |^'(c)|2sin2(c/2)dc, 
Jo 



inc/2 

and an orthonormal basis is the usual Fourier one. In) = -r=S= — — - with n G Z. 

' ' v47r sin(c/2) 

In this case the operator P = — ^(^ + 2 cot(c/2)), which eigenfunctions are 
|n), has eigenvalues Pn = ^n, the volume operator has eigenvalues 



^(n-l)/2 



and the operator c is well defined. In this case the eigenvalues of the inverse volume 
operator are difficult to calculate, but they can be approximated by El 

R / \ 15/2 

-^r;iiV2 27 when n « 2j 



dTn={ T^''{n-l)/2\2j) ^ ^ (Ig) 

>-l)/2 



^fniiV2 when n>2i, 



where j G N/2 is an ambiguity parameter. 

Now to obtain the effective Friedmann equation we use Bojowald's ideas (see 
for instance, [2 , 22 1): We identify the eigenvalues of p with a? to obtain the relation 

= ^n, then the eigenvalues of the inverse volume operator, as functions of a, 
can be written as follows: 

126 /' 3 \^^/^^12 „,u„„ „2 



d,(a)- ^Uj " ^hen a^«^j/?, ^^^^ 
y a ^ when a? » 7i/3. 

Then we change the Friedmann equation 
by the effective Friedmann equation 

■2 1 / 1 \ -3(l+a;) 

3 \aoJ ^ 

This equation coincides with the classical Friedmann one for large values of a, 
and thus, the big rip singularity (uj < —1) is not avoided because it occurs when 

a — )• 00. 
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On the other hand, when (oj > —1) there is a big bang singularity, which occurs 
when a — 0. In this case, near the singularity, the modified Friedmann equation 
models the classical one, but with the continuity equations P = ijJeffP where now 
We// = —5— 4a; < —1 because a; > —1, which means that the big bang singularity 
is avoided. 

To finish, note that all these calculation can be done in an easier way us- 
ing the /3 -representation assuming that (3 takes values in a circle. For example, 
choosing as a Hilbert space jC'^{—A, A), the operator V is given by — and 

then, the states |n) = '^'y^^ are orthonormal eigenf unctions with eigenvalues 
^j^. The parameter A could be determined identifying ^ with A^, which gives 
A = ^ = ^^3/4 - The eigenvalues of the volume and inverse volume operators 

are respectively Vn = X^\n\ and dn = + 2A|V2 _ _ 2A|V2||2^ 

then, performing the identification = X^n one gets 

. ^ / J^a^ when « 7A/2 
"^^""^ = [ when » 7A/2, ^''^ 

and now the modified Friedmann equation, when a; > —1, near the singularity 
models the classical one with the continuity equations P = ujeffP where ojgff = 
—3 — 2u! < —1, and thus, the big bang singularity is avoided, but not the big rip 
one. 

Remark 0.4. Other boundary conditions are possible in bounded domains, for 
example, homogeneous Dirichlet boundary conditions (physically it means that 
the Hubble parameter, or the connection, are into an infinite square well, which of 
course, does not have any physical basis). In this case, zero is not a eigenvalue of 
the volume, which allows us to define, without any ambiguity, the eigenvalues of the 
inverse volume operator as the inverse of the eigenvalues of the volume operator 

The problem in this formulation is that one has assumed that the values of the 
connection, or the Hubble parameter, belong in a circle rather in the real line (recall 
that in bounded domains a well-defined differential operator needs boundary con- 
ditions, in this case they have been assumed periodic, which means that the variable 
takes values in a circle). This assumption has no physical basis and for this reason 
it was left aside and it emerged the idea that the Besicovitch space could be the new 
Hilbert space in LQC. However as we have seen throughout the paper, this Hilbert 
space requires a new definition of the Hamiltonian operator, and the one used in 
LQC doesn't seem to be the correct one, because his classical analogue, after Len- 
gendre's transformation, gives rise to a non covariant Lagrangian. Moreover, using 
this Hamiltonian one obtains a modified Friedmann equations where the Hubble 
parameter is once again constrained to be bounded, and what is worse, it doesn't 
appear any higher-curvature term, like derivatives of the Hubble parameter, which 
are essential in modified or semiclassical cosmology. 
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A final remark is in order: In Quantum Gravity one usually follows the Dirac 
quantization scheme, that applied to cosmology leads to the well-known Wheeler- 
DeWitt equation. Appart from the unsolved problem of time, in LQC for the grav- 
itational part of the Hamiltonian is used the quantum version of ( 10 1, which as you 
have seen troughout the paper gives rise, at classical level, to a non covariant Lag- 
rangian. Then, all the conclusion obtained from this equation, particulary bounces 
(see for example |23]), need a deep revision. 



5. - Conclusions — Our main conclusion is that in loop quantum cosmology, con- 
trary to the current belief, the big rip singularity is not replaced by a non-singular 
bounce. What really happens is the theory is built from a Hamiltonian that, after 
a Legendre's transformation, gives a non covariant Lagrangian which is in contra- 
diction with the principles of General Gravity. Then, all the conclusions obtained 
from this Hamiltonian, in particular the big rip singularity, will be wrong. Another 
way to understand that the big rip singularity must survive in loop quantum cosmo- 
logy is the fact that this theory is based in the discrete nature of the space, which 
could influence the dynamics of our universe at small distances (at the order of 
Planck scale), but never at large distances where the discrete space could be subti- 
tuted by a continuous one, and thus, near the big rip the dynamics of the universe 
is unaffected by loop quantum corrections. 
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